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Where do the tedious products of C s come from? 

D. J. Broadhurst ^ 

^'Department of Physics and Astronomy, 

The Open University, Milton Keynes MK7 6AA, U.K. 

Lamentably, the full analytical content of the e-expansion of the master two-loop two-point function, with 
arbitrary self-energy insertions in 4 — 2e dimensions, is still unknown. Here we show that multiple zeta values 
(MZVs) of weights up to 12 suffice through 0(e®). Products of primitive MZVs are generated by a processes of 
“pseudo-exponentiation” whose combinatorics faithfully accord with expectations based on Kreimer’s modified 
shuffle product and on the Drinfeld-Deligne conjecture. The existence of such a mechanism, relating thousands 
of complicated rational numbers, enables us to identify precise and simple combinations of MZVs specific to 


quantum field theories in even numbers of spacetime 

1. Master two loop two point function 

The object of study in this contribution is the 
two-loop integral 

2{D/2-ae) r r 

I = j j ( 1 ) 

with og = 3Z?/2 —oi —02 —ct 3 —m —cts determin¬ 
ing the dependence on the number, D = A — 2e, 
of spacetime dimensions. It is clearly indepen¬ 
dent of the norm, p^, of the external momentum 
p. Remarkably, it has a 1440-element symmetry 
group [ |l| , corresponding to all permutations of 
6 linear combinations [ of the indices aj , com¬ 
bined with the total reflection aj D12 — aj. 

With parameters = 1 -I- rij-e, resulting from 
insertions on internal lines, the Taylor expansion 
of / = 6C(3) -I- 0(e) is known [ || to involve mul¬ 
tiple zeta values. 

The MZV C(5,3) = Em>n>o 
weight 8 and depth 2, appears in I at order 
and is also present in the six-loop beta function 
of (j)'^ theory [ H] . 

Similarly, C(3,5,3) = Ei>m>n>o 
with weight 11 and depth 3, appears at order e® 
and in cj)^ theory at 7 loops. 

In general, we expect L-loop counterterms to 
include MZVs of weights w < 2L — oi which 
the expansion of / through order is to be 

taken as a strong diagnostic. 


dimensions. 


It is one of the many scandals of our limited 
understanding of the analytical content of pertur¬ 
bative quantum field theory (QFT) that, despite 
many years of intense effort, we still do not know 
whether MZVs suffice for even the Taylor expan¬ 
sion of the two-loop integral 

Here we exhibit remarkable structure in the de¬ 
pendence of the X]fc=o Taylor coef¬ 

ficients through 0(e®) on the 

l + l + 2 + 2 + Z + A + h + l + ^ + l2 = A& (3) 

independent MZV structures from weight 3 to 
12. In particular, MZV structures that are im- 
primitive products will be generated by a process 
of “pseudo-exponentiation” whose very existence 
appears to me to be little short of miraculous and 
whose details resonate strongly both with Dirk 
Kreimer’s combinatorically modified shuffle prod¬ 
uct, discussed in these proceedings, and also with 
the Drinfeld-Deligne conjecture, namely that the 
Grothendieck-Teichmiiller algebra has precisely 
one generator of each odd degree greater than 
unity, and none of even degree. 

2. Choices of primitive MZVs 

As Don Zagier observed, the Drinfeld-Deligne 
conjecture leads to an enumeration P{n) = P(n — 
2) -b P{n — 3) of Q-linearly (i.e. rationally) inde¬ 
pendent MZV structures of weight n, seeded by 
P(l) = 0, P{2) = P(3) = 1. This generates the 
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Padovan numbers P(3) to -P(12) in (||). The fil¬ 
tration of products into primitives then requires 

M(w) =-'^n{w/n)Q{n) (4) 

w ^ 

'n\w 

primitives to be chosen at weight w, where 
Q(n) = Q{n — 2) -I- Q(n — 3) is the Perrin se¬ 
quence, seeded by <5(1) = 0, Q(2) = 2, Q(3) = 3, 
and the Mobius function ^ ensures that 

i(5) 
I — — X'^ 

n>0 

At weight w = 2 we have the primitive tt^ (which 
mathematicians sometimes relegate to re = 0) 
and at each odd weight 2n -I- 1 > 1 we have the 
primitive (^{2n + l). That leaves us with one more 
primitive MZV to choose at each of the weights 8, 
10 and 11. We might take these as ^(6, 2), ^(8, 2) 
and C(8, 2,1). At weight ic = 12 we need a primi¬ 
tive of depth 2, such as C(10, 2), and one of depth 
4, such as C(4,4, 2, 2). 

From a mathematical point of view, these are 
entirely adequate choices of primitives, since ev¬ 
ery MZV of weight re < 12 is expressible as a 
rational linear combination of products of the fil¬ 
tration. 

However, we are concerned here with QFT, 
where more specific combinations of MZVs are 
found to occur in counterterms [ ^ . I shall sketch 
how I was led, by Kreimer’s fertile ideas [ 
and by the Drinfeld-Deligne conjecture, to the 
following set of QFT primitives up to weight 12. 

1. At weights w = 2,3, 5, 7,9,11, take ({w). 

2. At w = 8, take ^(5,3). 

3. Atw;=10, take C(7,3)-k3C^(5). 

4. Atw = 11, take C(3,5,3) - C(3)C(5,3). 

5. Atw= 12, take ({7, 5) and the combination 
of alternating Euler sums below. 

The choice of a second primitive at weight 12 is 

Z2(12) = C(7,5) - C(5,7) + C(5)C(7) (6) 


comprising Euler sums of the form 

= E (7) 

^ ^ m°' n" 

m>n>0 

C(«) = E ^ (8) 

^' n 

n>0 

with alternation of sign notated by a bar above 
the corresponding index. 

With these very specific choices, products in 
the Taylor expansion of (|^) through order e® 
will be generated by pseudo-exponentiation. The 
choice of primitives which makes this possible is 
essentially unique. The existence of such a choice 
is truly remarkable. For example, the choice of 
primitive in (i was determined by a tiny subset 
of the available data, whereupon transformation 
to the appropriate basis reduced a file of many 
megabytes of tedious products of C’s to a single 
line of code. 


3. The master function 


The key to the expansion of (|^) is the function [ 

I 


S{a,b,c,d) = 


TT cot TTC 


1 


Hla,b,c,d) c 

^ ^ F{a + c, -b, -c, b + d) (9) 
be 


with a hypergeometric series of the form 

( 6)77 


F{a,b,c,d) = ^ 


(1 -f 0)71(1 -|- d)r 


( 10 ) 


where {a)n = r(a -I- n)/r(a), and with a ratio of 
products of Gamma functions of the form 


H{a, b, c, d) 


G{{a,b,c,d,a + b + c + d}) 
G(^\^a c, a d, b A c, b d}) 


where G{S) = OaeS ^(1 + a). 

If the propagators in two adjacent internal lines 
of the master two-loop diagram have unit expo¬ 
nents (i.e. no dressing by self-energy insertions) 
then the recurrence relations from integration by 
parts may be solved in terms of two instances [ ^ 
of the construct (||), which solves 


aS{a, b,c,d) = l + 


(a + c) (a + d) 


S{a 


1,5, c, d) 


a A b A c A d 
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in the presence of the symmetries 

S{a, 5, c, d) = S{b, a, c, d) = —S{c, d, a, b). (12) 

Moreover, the 1440-fold symmetry of the mas¬ 
ter diagram enables one to reconstruct the 9th- 
order Taylor expansion of /, with 6 arguments, 
from the llth-order expansion of 5", with only 4 
arguments. If there is any new type of number, 
beyond the MZVs that are generated Iw the e- 
expansions of Pochhammer symbols in (|l0| ) , then 
it will show up only at weights greater than 12. 

4. Pushdown at weight 12 

The technology for developing the Taylor ex¬ 
pansion of S through weight 11 was developed 
in [ H. At weight 12 we now encounter a remark¬ 
able phenomenon: the depth-4 MZV C(4,4, 2,2) 
is primitive within the confines of the theory of 
MZVs, but widening the analysis to include al¬ 
ternating Euler sums, one discovers that it may 
be eliminated in favour of a depth-2 sum, such as 
C(9,3) = Specifically, 

the combination C(4,4, 2,2) — (8/3)^C(9 j 3) is re¬ 
ducible to 11 weight-12 terms in the MZV basis, 
thus showing that C(9, 3) can furnish the 12th. 

Here one sees another scandal of our limited 
understanding. This “pushdown” from depth- 
4 MZVs to depth-2 Euler sums was found em¬ 
pirically using David Bailey’s implementation 
of Helaman Ferguson’s PSLQ algorithm, which 
turns high precision numerical evaluations into 
immensely probable integer relations. In this 
case, the effort was infinitesimal, in comparison 
with the PSLQ investigations in [ |^. Never¬ 
theless, no-one seems to have any idea of how 
to prove such pushdowns rigorously, apart from 
labouring to solve the huge complex of linear re¬ 
lations between 3^^ = 531441 words of 12 let¬ 
ters, taken from the alphabet A = dx/x, B = 
da:/(l — x), C = —dx/(l -f x) for the iterated 
integrals corresponding to Euler sums. 

In terms of CPUtime, it takes but a sec¬ 
ond of PSLQ work to find the simple coefficient 
(8/3)^ of C(9, 3) in the pushdown of ^(4,4, 2,2). 
Yet it would appear to require a huge expense 
of computer algebra, feasible perhaps only with 
Jos Vermaseren’s form, to obtain this reduc¬ 


tion by systematic application of the shuffle al¬ 
gebras of nested sums and iterated integrals in 
the {A,B,C} alphabet [ p| 

To develop the Taylor expansion of S to weight 
12, PSLQ was used to find 139 such reductions of 
MZVs to a basis involving C(9,3). 

5. Pseudo exponentiation 

After achieving the pushdown to depth 2, one 
is left with something that seems hugely indi¬ 
gestible, namely the reduction of (^ 3 ^) = 

1364 Taylor coefficients of S to Q-linear combi¬ 
nations of P{k -f 1) = 47 terms formed by 
products of primitives. The data comprises many 
thousands of rational numbers which are typically 
the ratios of 10-digit integers. Some reduction 
of this complexity is achieved by using the sym¬ 
metries (p^. For example the leading term in 
S{a, b, c, d) is simply (a -I- 6 — c — d)({2) and this 
cancels when one combines the two instances of S 
in I — 6^(3) -I- 0(e). Even after taking account of 
such obvious symmetries, several thousand prod¬ 
uct terms remain to be explained. 

Discussion with Dirk Kreimer suggested that 
perhaps some of the tedious products of primi¬ 
tives might follow a pattern. The basic idea was 
that some simple structure might be delicately 
“composed” with an expression relatively free of 
such products, so that, for example, one might 
write S = E V R, where E is made of Gamma 
functions (and hence exponentiates a series of 
^’s), R is far simpler than S, and the rule of com¬ 
position, denoted by V, has the distinctive QFT 
property of counting like objects and then modi¬ 
fying the ordinary process of multiplication by an 
appropriate symmetry factor. 

There was an obvious candidate for E: the re¬ 
ciprocal of E[ in the first term of (^, correspond¬ 
ing to exponentiation of 

-logJ7(a,6, c,d) = 

, s 

s>l 

N{s) =a^+E Ed'^ + {a + b + c + dy, 

D{s) = (a -I- c)® -I- (a -I- d)® -I- (6 -I- c)® -I- (6 -I- d)®. 

However it soon became clear that the powers 
of TT^ in log id can play no role in the composi- 
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tion Ey R. This seemed to chime well with the 
Drinfeld-Deligne conjecture. Hence the Ansatz 


E{a, 5, c, d) 


V H{a,b,c,d) ) 


(13) 


was made, since it exponentiates only odd C’s. 

The challenge was then to find a product-free 
expansion of the form 


12 

R{a,b,c,d) = E C{w)Jnj{a,b,c, d) 

w^2 

+ ^ Z{w)Kuj{a,b,c,d) 

w=8,10,ll 

+ Zk{l2)Ki2Aa,b,c,d) (14) 

fc=l,2 


involving the 11 single sums ({w) = X)n>o 
with w G [2,12], a single additional primitive 
MZV at each of the weights w = 8 , 10 , 11 , and 
two such primitives a,t w = 12 . Since all 16 of 
these terms are present in S, and E = 1 + 0{e^) 
does not exponentiate powers of tt^, this is the 
minimal product-free Ansatz for R. What was 
far from clear is that there exists a well-defined 
procedure Ey R that can generate the remaining 
47—16 = 31 product structures in S from the 
11 + 5 = 16 single-sum and primitive terms in R. 
Bearing in mind that a weight-w product term in 
S entails all monomials of degree w — 1 in 4 

variables, the odds against finding a composition 
procedure appeared at first so overwhelming as 
to make such an enterprise foolish. 

However, following the well tried maxim that 
“QFT is far smarter than we are”, I looked to see 
what happens when one takes the composition 
to be that of simple multiplication. This faith¬ 
fully reproduced product terms such as ^( 3 )^( 2 ) 
and C(3)C(4), but to Dirk Kreimer’s great delight 
gave the wrong result for [^( 3 )]^, whose true de¬ 
pendence on the 56 monomials of 5th order in 4 
variables is proportional to (A^(3) — D{3))J3 yet 
has precisely half the coefficient that would be 
predicted by straightforward multiplication. In 
other words the composition E y R must take 
account of symmetry factors: if [('(n)]^/p! from 
the exponential E is composed with C,{n) from 
the product-free expansion R, the contribution to 


Ey R must be taken as [C,{n)Y^^ / {p -I- 1)!. Thus 
to generate the 165 monomials in the [C(3)]^ term 
of S we must modify the result of naive multipli¬ 
cation by the factor 2!/3! = 1/3. That aspect of 
pseudo-exponentiation had been eagerly antici¬ 
pated. 

A second aspect came as a surprise, though in 
retrospect it signals that QFT is in deep accord 
with the Drinfeld-Deligne conjecture. The first 
example occurs at weight 8, where we know that 
we must include a new primitive, such as ('(5,3), 
in Z(8). To a mathematician, such a primitive 
could include an arbitrary rational multiple of 
C(5)C(3). But we seek to generate all monomi¬ 
als found to be multiplying C(5)C(3) in S, using 
the composition Ey R. There are, at most, three 
sources: C(5) V C(3), C(3) V C(5) and, possibly, a 
multiple of C(5)C(3) from 1 V ^(8). The func¬ 
tional dependencies on a, b, c, d of all three con¬ 
tributions are predetermined; we have just three 
rational numbers at our disposal to fit 120 mono¬ 
mials. A solution exists and it is unique: every¬ 
thing is given by ({3) V (/(5). 

It appears that the QFT expansion “knows 
more than we do” about the structure of MZVs. 
Faced with a choice between ^(5) V C(3) and 
C(3) V C(5) it rejects the former and takes naive 
multiplication as the composition rule for the lat¬ 
ter. For good measure it tells us that the new 
primitive is Z{8) = ^(5,3), modulo tt®, with no 
contribution from C(5)/'(3). Note that we would 
get an unambiguously wrong answer were we to 
take the primitive as C(3, 5). 

The situation becomes yet more interesting 
at weight 10, where ({7) V f{3) is rejected, 
C(3) V C(7) is accepted, as a plain multiplication, 
and the primitive is determined to be .^(10) = 
C(7, 3) -I- 3C^(5), modulo 7r^°, with no possibil¬ 
ity of /’(7)/’(3) appearing in this primitive. Here 
one sees the monomials in the [C(5)]^ term of S 
generated by the sum of two terms: {N{5) — 
D(5))J5 and Kio, each of which was predeter¬ 
mined. The former is normalized, as expected, 
by taking account of the factor of 1/2 in pseudo¬ 
exponentiation; for the latter we have only a sin¬ 
gle rational number, with which to fit 220 mono¬ 
mials, and the coefficient 3 in C(7, 3) -I- 3C^(5) 
emerges in every case. 



5 


At weight 11, the working out is even more 
impressive: [C(3)]^ V C(5) and ^(3) V C(5i3) are 
present, as naive multiplications; (C(5)C(3))VC(3) 
is rejected; and the primitive is given for each of 
the 286 monomials as the combination ^(11) = 
C(3, 5,3) —C(3)C(5,3), in precise accord with anal¬ 
ysis of subdivergence-free 7-loop counterterms [ 

I- 

6. Experimentum crucis 

Finally, QFT seems to “know more than we 
do” about pushdown of depth-4 MZVs to depth- 
2 alternating Euler sums at weight 12. There 
still remain several hundred rational numbers to 
be generated by a procedure that is now unam¬ 
biguous: [C(3)]^ VC(3) = [C(3)]^/4 is combinator- 
ically modified, a la Kreimer, and ((7)\/((5) = 0, 
C(9)VC(3) = 0, are rejected, a la Drinfeld-Deligne, 
in favour of a pair of primitives. Taking the first 
as Zi(12) = ({7, 5) we then determine the second, 
^ 2 ( 12 ), modulo the first, and modulo tt^^. It is 
forced to contain all 10 of the remaining MZV 
structures in a basis that includes C(4,4,2,2). 
Moreover the coefficients of these terms have nu¬ 
merators of typically 10 digits. 

When we push C(4,4, 2,2) down to depth 2, 
via its empirical PSLQ relation to C(9, 3), some 
sanity emerges, since the terms in .^ 2 ( 12 ) now in¬ 
volve only double sums and products of two sin¬ 
gle sums. Yet still there are apparently grotesque 
rational coefficients. Finally, we discovered that 
this seemingly random assortment of terms is 
simply given by (|^), with 3 terms whose coeffi¬ 
cients are ±1, though the simplicity of C(7)C(5) = 
945C(7)C(5)/1024 would have been disguised had 
we not transformed to a product of alternating 
sums. 

It is as if QFT were taunting us with our ig¬ 
norance of the mapping between diagrams and 
numbers that results from the Feynman rules. A 
vast quantity of data, collected by painfully inad¬ 
equate methods, collapses to an amazingly sim¬ 
ple answer. We are, physicists and mathemati¬ 
cians alike, stumbling on the edge of a structure 
that is far more refined than the clumsy meth¬ 
ods by which we investigate it. In this partic¬ 
ular instance, it was the sheer volume of appar¬ 


ently random - yet in fact gracefully proportioned 
- data, from many numerical instances of self¬ 
energy dressings, that led to an answer to the 
question posed in the title of this contribution. 
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